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Abstract 

We examine decays of a spin-1 bottomonium into a pair of light scalar Dark Matter (DM) 
particles, assuming that Dark Matter is produced due to exchange of heavy degrees of freedom. We 
perform a model-independent analysis and derive formulae for the branching ratios of these decays. 
We confront our calculation results with the experimental data. We show that the considered 
branching ratios are within the reach of the present BaBaR experimental sensitivity. Thus, Dark 
Matter production in T decays leads to constraints on parameters of various models containing a 
light spin-0 DM particle. We illustrate this for the models with a "WIMP less miracle", in particular 
for a Gauge Mediated SUSY breaking scenario, with a spin-0 DM particle in the hidden sector. 
Another example considered is the type II 2IIDM with a scalar DM particle. 
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1. INTRODUCTION 



The origin, presence and nature of Dark Matter (DM) in our Universe remains one of 
the biggest mysteries of the particle physics and astrophysics [jj. Understanding the nature 
of Dark Matter, which accounts for majority of mass in the Universe, represents a crucial 
step in connecting astronomical observations with predictions of various elementary particle 
theories. 

Many theories, the extensions of the Standard Model (SM), predict one or more stable, 
electrically-neutral particle (s) in their spectrum, which can be possible Dark Matter candi- 
date(s). Different models provide different assignments for DM particles' spin and various 
windows for their masses and couplings. To test this great variety of hypotheses, several 
techniques for DM direct or indirect search are currently developed. 

Recent experimental measurements of Dark Matter relic abundance, Qj^Mh"^ ~ 0.11 by 
WMAP collaboration ^, can be used to place constraints on the masses and interaction 
strengths of DM particles. Indeed, the relation 
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with M and g being the mass and the interaction strength associated with DM annihilation, 
implies that, for a weakly-interacting particle, the mass scale should be set around the 
electroweak scale. This, coupled with an observation that very light DM particles might 
overdose the Universe (known as a Lee- Weinberg limit 3|) , seems to exclude the possibility 
of light weakly- interacting massive particle (WIMP) solution for DM, setting Mj^m > 2 — 
6 GeV. 

A detailed look at this argument reveals that those constraints can be easily avoided in 
concrete models, so even MeV-scale particles can be good candidates for DM. For instance, 
low-energy resonances - such a light CP-odd Higgs in the MSSM extensions with a Higgs 
singlet, or a light extra gauge U-boson - could enhance the DM annihilation cross-section, 
without the need for a large coupling constant [sl, jfil-jsi. Even if no light resonances 
exist, the usual suppression of DM annihilation cross-section, M^^/M^, used in setting the 
ee- Weinberg limit, does not hold, if Dark Matter is a non-fermionic (e.g. spin-0) state 
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151] ■ Furthermore, DM production could be non-thermal 16|, llTI], in 



which case the constraint provided by Eq. fillip does not apply. Thus, models with light 
Dark Matter {Mdm ~ few GeV or less) still deserve a detailed and thorough study. 



In this paper we consider the possibihty of using of T meson decays with missing energy, 
to test the models with a hght spin-0 DM particle. Studies of heavy meson (and in particular 
T meson) decays with missing energy may be especially valuable in light of the fact that 
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DM direct search experiments, such as DAMA 
which rely on the measurement of kinematic recoil of nuclei in DM interactions, lose (for cold 
DM particles) sensitivity with decreasing mass of the WIMP, as the recoil energy becomes 
small. Of course, light Dark Matter may also be produced at high energy colliders, however 
the production rate is naturally going to be insensitive to precise value of the WIMP mass, 
if the one is much less than the beam center-of-mass energy 2J|. Indirect experiments, 
such as HESS [25], are specifically tuned to see large energy secondaries, only possible 
'or weak-scale WIMPs._The backgrounds for positron and antiproton searches by HEAT 

281 ] experiments could be prohibitively large at small energies. 
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Thus, T (and/or other heavy meson) decays with missing energy may serve as alternative 
DM search channels, capable to provide us with an information on the WIMP mass range, 
hardly testable by the experiments discussed above. 

So far, T meson decays into Dark Matter have been considered within the models, where 
DM partic 
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es interaction with an ordinary matter is mediated by some light degree of free- 
2^. Apart from desire of having DM annihilation enhancement (due to a light 
intermediate resonance) and thus having no tension with the DM relic abundance condi- 
tion, it is also known that T meson SM decay is predominantly due to strong interactions. 
Thus, the WIMP production branching ratio, in general, is greatly suppressed compared to 
relevant weak B decays, and in particular to B —>■ K + invisible transition UJ, Il2|]. In light 
of this, it might seem natural to concentrate only on the models within which Dark Matter 
production in T decays is enhanced due to exchange of a light particle propagator. 

Yet, the aim of the present paper is to study T(IS') decay into a pair of spin-0 DM 
particles, T(IS') — > $$*, and T{3S) decay into a pair of spin-0 DM particles and a photon, 
T(3S') $$*7, within the models where light Dark Matter interaction with an ordinary 
matter is due to exchange of heavy particles (with masses exceeding the bottomonium mass). 
As it is mentioned above, these models may be free of tension related to satisfying the DM 
relic abundance constraint as well. Examples of such models will be discussed below. Also, 
new experimental data on T decays into invisible states have been reported by the BaBaR 
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collaboration [30|, l31|- According to these data, 

B{T{1S) invisible) < 3 x 10"^ ([112) 

and 

5(T(3S) 7 + invisible) < (0.7 - 31) x 10"*^ (^3) 

where the interval in the r.h.s. of Eq. fll[3p is related to the choice of the final state missing 
mass. These bounds are significantly stronger than those on invisible Tfl^*) decays (with 
or without a photon emission), reported previously by Belle and CLEO 
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by Particle Data Group [3^. We show here that BaBaR experimental data on T meson 
invisible decays (without or with a photon emission) may constrain the parameter space of 
light scalar Dark Matter models, even if there is no Dark Matter production enhancement 
due to light intermediate states. 

We also illustrate (in Sections [5] and [6]) that the study of Dark Matter production in 
T decays allows us to test regions of parameter space of light spin-0 DM models that are 
inaccessible for B meson decays with missing energy. It is also worth mentioning that T 
decays are sensitive to a wider range of WIMP mass than B decays. Thus, the study of 
WIMP production in T decays is complementary to that for B meson decays. 

Within the models where scalar DM consists of particles that are their own antiparticles, 
only T $$7 transition is relevant. The transition rate of T — $$* vanishes, if $ = $*. 
Indeed, by angular momentum conservation, the produced DM particles pair in T — > 
decay must be in a P-wave state, which is impossible because of the Bose-Einstein symmetry 
of identical spin-0 particles. 

The scenarios with light complex scalar Dark Matter, albeit implying some continuous 
symmetry related to the internal charge of the complex (electrically) neutral scalar field, 
may be realized in many extensions of the Standard Model. Some of the models allow the 
scenarios both with a real and with a complex light scalar DM field {gI-IiS. 
these models, study of the decay T(15') $$* represents an excellent opportunity to test 
the DM field nature either at present or in the future (if higher experimental resolution is 
needed) . 

As mentioned above, we restrict ourselves to the class of models where light spin-0 DM 
production is mediated by heavy degrees of freedom. At the energy scales, associated with 
T decays, heavy intermediate degrees of freedom may be integrated out, thus leading to 



14 



15 



35|. Within 



a low-energy effective theory of four-particle interactions. Our strategy would be deriving 
first model-independent formulae for the T(15') —>■ $$* and T(3S') $$*7 branching 
ratios within the low-energy effective theory. Then, we confront our predictions with the 
experimental data, deriving model-independent bounds in terms of the Wilson operator 
expansion coefficients as the parameters that carry the information on an underlying New 
Physics (NP) model. Finally, within a given model, using the matching conditions for the 
Wilson coefficients, we translate these bounds into those on the relevant parameters of the 
considered model. 

The paper is organized as follows. In Section [2] the general formalism is developed and 
model-independent formulae for the T(IS') — > $$* and T(3S') — > $$*7 branching ratios are 
derived. Both the case of complex and the case of self-conjugate DM field are considered. In 
Section [3] the neutrino background effect is discussed. We show that it has negligible impact 
on our analysis. In Section H] we confront our calculation results with the experimental data 
and derive model-independent bounds on the Wilson coefficients. In the next sections we 
transform these bounds into constraints on the parameters of particular models. In Section [5] 
models with a complex spin-0 DM field are considered. We choose the class of mirror fermion 
models as an example. These models include, in particular. Gauge Mediated SUSY Breaking 
scenarios with the DM particle in the hidden sector and the mirror fermions being connectors 
between the hidden and the MSSM sectors [l3, Q, [l^. An example of the self-conjugate 
DM scenario, two-Higgs doublet model with a scalar DM particle, is considered in Section [6l 
The concluding remarks are given in Section [71 

2. GENERAL FORMALISM: MODEL-INDEPENDENT FORMULAE FOR THE 
DECAYS BRANCHING RATIOS 

We treat T states - neglecting the sea quark and gluon distributions - as bound states 
of bb valence quark-antiquark pair that annihilates - with or without emission of a photon - 
into a pair of Dark Matter particles. To this approximation, the relevant low-energy effective 
Hamiltonian may be written as 

Heff = ^T.C^O, Pl) 



where Ah is the heavy mass and 
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Oi = [b b) ($*$) , O2 = tmi, (6756) ($*$) , 

O3 = (hrb) , O4 = (&7^75&) ($*^^^$) P2) 

with 9 = l/2{d — d). It is worth noting that with the notations used in (121 ip and 
(121 2p . all the operators Oi, i=l,..4, are Hermitean, thus all the Wilson coefficients Ci must 
be real. Another point to be made is that the Ci are low-energy renormalization scale- 
independent. This stems from the renormalization scale invariance of the hadronic parts of 
operators Oi, which - combined with the fact that DM particles do not interact strongly or 
electromagnetically - leads to low-energy scale invariance of Oi and (from the scale invariance 
of the effective Hamiltonian) to that of Ci. If DM consists of particles that are their own 
antiparticles, then only first two operators in Eq. f l2[2p would contribute. 

To the considered approximation, the operator basis, given by Eq. fl2[2p . is the most gen- 
eral one for the T(IS') $$* transition. With the use of the same basis, the other decay, 
T{3S) — > $$*7, occurs by means of bi- local interactions depicted in Fig. [1] (a), (b). In 
principle, it is also possible - if the intermediate heavy states are (electrically) charged - that 
both the photon and the DM particles originate from the same effective vertex (Fig. [11(c)). 
In order to take into account this diagram, higher dimension operators, involving electro- 
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magnetic field tensor, should also be included. However, such higher dimension operators 
would enter the expression for H^ff with higher powers of 1/Ah factor. Extra suppression 
of diagram in Fig. [T](c) in powers of the heavy mass inverse may be seen from the fact that 
one must introduce an extra heavy propagator or, equivalently, an extra power of 1/Ah or 
1/A^, as the photon is emitted by a heavy intermediate degree of freedom. 

Thus, one may neglect contribution of the diagram in Fig. W^c), as compared to that of 
the other two diagrams. This justifies the use of the operators basis, given by Eq. fl2[2l) . for 
T{3S) $$*7 decay as well. 

In the rest frame of T meson, matrix elements of the local hadronic currents may be 
parameterized at the origin, x = 0, as 

(0| 6(0)6(0) |T) = (0| 6(0)756(0) |T) = (0| 6(0)7^75^(0) |T) = 

(0| 6(0)7^6(0) |T) = frMrei^ip), (0| 6(0)a'^^6(0) |T) = -tfr [ p", 4(p) ] (03) 

where /x is the decay constant, Mx is the mass, p = (Mx,0) is the momentum and ex(p) 
is the polarization vector of T meson. Although the tensor current bcr'^'^b is not present 
in fl2[2l) . it inevitably appears in calculations of the T{3S) $$*7 amplitude. Thus, we 
need this current matrix element parametrization as well. Also, it is worth mentioning that 
p ■ er = due to the vector current conservation. 

Non-local hadronic currents matrix elements may be expressed, to the leading order in 
l/mh expansion, in terms of those of the local currents in the following way: 

(0| 6(xi) r 6(x2) |T) = e-^(p/2) (-i+-2) 6(0) F 6(0) |T) (04) 

where F is a product of 7-matrices. This relationship is derived using the constituent quark 
approach, which assumes that within T meson both 6 and 6 are static and have a mass 
Mr/2 each. Thus, one neglects 0{Aqcd) difference between Mr/2 and rrib as well as the 
quark- ant iquark Fermi motion effects. 

Let us proceed to the transition amplitudes and rates. As it is discussed above, to the 
leading order the decay T(IS') $$* occurs as a result of an effective local four-particle 
interaction. Thus, as it follows from fl2[ll) . fl2[2l) . fl2[3l) . only operator O3 contributes to this 
decay. Furthermore, if $ is a real scalar (or pseudoscalar) state, it is easy to show that 
contribution of O3 vanishes as well. 
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Calculation of T{1S) $$* branching ratio is straightforward: for $ being a complex 
scalar state one gets 



r(T(l^) ^ CI /■ 



3/2 



J- T(15) A|f 487ri T{15) L ^ ^ J 

where m$ is the DM particle mass and rx(is) = (54.02 ± 1.25)keV js^ is the T(15') total 
width. 

For $ being a self-conjugate spin-0 state, $ = $*, one has 

5(T(15) ^ $$) = (06) 

As it was mentioned above, this result is related to the fact that the final DM particle pair 
state must be a P-wave, which is impossible due to the Bose-Einstein symmetry of identical 
spin-0 particles. In what follows, r(T(15') —>■ $$) must also vanish in higher orders in l/rrib 
operator product expansion. 

Thus, provided that DM pair production is the dominant invisible channel, the signal 
for T(IS') — s> invisible decay would imply that the light spin-0 DM field has a complex 
nature. No evidence for the T(IS') invisible mode allows one to put some constraints 
on the parameters of the models with light complex scalar Dark Matter, as we illustrate in 
Sections H] and [5l 

Consider the other mode, T(3S') —>■ $$*7. As it was mentioned above, this decay occurs 
as result of a bi-local interaction, and as direct calculations show the contribution of operator 
O3 to the decay amplitude is vanishing, whereas the other operators have, in general, a non- 
zero contribution. In the case of $ being a complex scalar state, one gets 

A{T{3S) $<l>*7) = Ai{T{3S) <l>$*7) + A2(T(35) ^ $<l>*7) + 

+A^{T{3S) $<l>*7) (07) 



where 



Ai(T(35) ^ <l>$*7) = ^^^^ e^ri3S){p) [Mujg,.-p,K] (08) 



A2(T(35) - $$*7) = ^^^^ ^M-A P6*"(fc) fe\^,s){p) P9) 
A,{T{3S) -> <l>$*7) = -^^^^ ^^M-A Fe*°(A;) (p2 - Pi)^e^(35)(p) PlO) 
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Here k = {u, k) is the photon momentum, e{k) is the photon polarization vector, pi and p2 are 
the momenta of the DM particle and antiparticle respectively (by momentum conservation, 
Pi + P2 + k = p). 

Note that there is a 5 = 7r/2 difference in the phases of A2 and A^. Thus, these parts 
of T(3S') $$*7 amplitude do not interfere. It is also easy to check, after doing some 
algebra, that there is no interference between Ai and A2 or A^ as well. This is also easy to 
understand: contribution of the parity-conserving operator Oi does not interfere with that 
of the parity violating operators O2 and O4. 

In what follows, the differential branching ratio for T(3S') $$*7 decay may be written 
in the following form: 

^(T(35) - <|.$*7) = Tr^^(T(3^) - $$*7) = 
ds rT(35') ds 

JTD JTD JJD 

= -j^{T{?>S) ^ $<l>*7) + -#(T(3^) ^ $$*7) + -77^(T(3^) ^ <l>$*7) (011) 
as as as 



where rT(35) = (20.32 ± 1.85)fceF y is the T{3S) total width and 



^(T(3^) - $$*7) = (012) 
ds ^ ^ ' ' A.\i 47r 547rrT(35) V s ^ ^ 

^(T(3S) - $$*7) = Trf ^ ^ (013) 

ds A|^47r 1627rrT(3s) \ ^ ) 

Here a = 1/137 is the electromagnetic coupling constant, x$ = '^|/-^t(35) and s = 
s/M|-(3_5-), where 

s = (Pi + P2)' = (p - A;)' = M2(3^) - 2 c^AfT(35) 

is the missing mass squared. Note that the kinematically allowed range for the missing mass 
squared is 4m| < s < My(35). Subsequently, 4a;<j> < s < 1. 

The partially mie^rated branching ratio for T(3S') — > $$*7 decay. 



S(T(35) ^ <l>$*7)|.<s_ = -r(T(35) ^ $$*7) Pl4) 

J ix^ as 

with Smaa; < 1 or, equivalcutly, Smax < ^T(35)' is given by 

i?(T(3S) ^ $<l>*7)|s<._ = i?i(T(35) ^ $$*7)|s<._ + 

52(T(35) ^ *7)|s<._ + i?4(T(35) ^ $ $*7)|.<s_ Pl5) 
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where 



5i,2(T(35) ^ <l><l>*7) 



I S'^^max 



C?2 a 



P 



T{35) 



Ajj Atx 1087rrT(3S)Mx(35) . 

,2 



2M- 



T(3S) 



+2mi 



-Ami) - 



T{35) 



In 



+ 



4m| 



i?4(T(35) ^ $$*7)|.<.„ 



T(3S) 



K\j All 1627rrT(3S)^T{35) 
' ^m.n.x \ 9 



2md 



(016) 



3 A/f2 



3^T(3S) , 
't(35) 
^'^T{35) 



2ml 

/lf2 

^^^T{35) J 



max \ '^max 



Amt 



Smax + \ S. 



Ami 



2ma 



(017) 



We use the partially integrated branching ratio to confront the theoretical predictions with 
experimental data. The existing experimental bounds on T — 7 + invisible mode are 



final state missing mass dependent 3ll, 133|. In particular, they may be very loose, or even 
there may be no bound, when the missing mass is close to its upper threshold, the T mass. 
Thus, imposing a cutoff s < Smax enables one to use existing experimental constraints on 
T 7 + invisible transition in the most efficient way. Also, this way one gets rid of the 
missing mass range, where the emitted photon energy is 0{Aqcd) and non-perturbative 
QCD effects may be of importance. The price we pay is shrinking of the WIMP mass range, 
where our analysis is efficient. 

Taking Smax = ^t(35)' S^^s total integrated branching ratio for T(3S') $$*7 
decay: 

B{T{3S) $$*7) = 5i(T(35) ^ <l><l>*7) + 52(T(3S) -> <l><l>*7) + 

+54(T(3^) ^ <l>$*7) 

where 

'■'^ ^ ^ Af, 47r 1087rrx(35) 



P18) 



(^T(35) "^1 y 



(T(3S) ^ = 



fr{3S) 



Mt(3S) 

2 31 2 3m 



T(3S)^ 





-Aml- 




-4m|Y 


2m$ j 




-Ami - 



(019) 
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6m| (3M^(3^) + 2m| 



In 



2m<j, 



4m| 



(020) 



The total integrated branching ratio may be used for theoretical studies, in particular to 
specify if a given model or class of models may be tested with the existing level of experi- 
mental accuracy. 

In the case of $ being a real scalar (or pseudoscalar) state, contribution of O4 vanishes 
and contributions of Oi and O2 to the branching ratio must be multiplied by factor two. In 
this case, one can rewrite formulae fl2[lll) - fl2[20p in a more simple form: 



dB 
ds 



(T(35) ^ $$7) 



An 



27'kT 



T(35) 



(021) 



i?(T(35) ^ $$7)|.<.„ 



{CI + CI) a 



T(3S) 



A% An 547rr 



T{35)^T(3S) 



+ 



+2m| 



-Ami] - 



+ \ s 



- Ami 



2ma 



P22) 



and 



5(T(35) ^ $$7) 



{CI + C|) g /|(35) 

47r 547rrx(35) . 



T(3S) 



+ 2m: 



8m| (M- 



T(3S) 



In 



T(3S) 



Mx(35) + V^M|(3g) - Am% 
2m* 



(023) 



As it was mentioned above, T(3S') ^ $$7 is the only mode that we use to test the models 
with self-conjugate spin-0 light Dark Matter. 



3. NEUTRINO BACKGROUND 



In this section we consider neutrino-antineutrino pair production within T(IS') — *■ 
invisible and T(3S') — ^ 7 + invisible channels. We examine possible impact of the neu- 
trino background on our analysis. 

Neutrino background in T(IS') invisible and T(3S') — ^ 7 + invisible decays occurs 
due io T {IS) ^ vv and T(3S') uu'-f transitions. These transitions may occur both within 
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the Standard Model and due to New Physics interactions. The NP contribution is model 
dependent and should be examined (upon necessity) in the framework of a particular model. 
Here we concentrate on the Standard Model contribution only. To the leading order within 
the SM, bb —>■ vv transition is mediated by virtual Z boson. 

T ^ z/i/ and T — * vv^ transitions have been originally discussed in [36] . Later on T — > z/i/ 
has been studied in detail in js^] , both within the Standard Model and within of some of its 
extensions. The SM result of ref. jsT] may be rewritten in terms of T(IS') vv branching 



ratio as 

B{T{1S) ^v-v) = ^'^'-> " " ^ = ^ ( 1 - sin^ Qw ) ' '''l^" '^'^^ pi) 

i T(is) 4»7r \ d / i x{i5) 

where Gp vs, the Fermi coupling, Qy^ is the weak mixing angle, and = 3 is the number 
of light non-sterile neutrinos. We use Gp = 1.166 x 10~^GeV~'^, sin^ 6w = 0.231, Mr(is) = 



9A5GeV and rT(is) = (54.02 ± l.25)keV 



34]. The decay constant /t(is) niay be extracted 



from the experimental measurements of T(IS') — e~^e rate: one gets /t{i5) = (0.715 ± 



0.005)GeV^ [381]. Using indicated values of the input parameters, one gets 

S(T(15) ^ z/z?) = (1.03 ±0.04) X 10"^ ([12) 

Thus, T(IS') vv decay branching ratio is about 30 times less than the BaBaR experi- 
mental bound on i?(T(lS') —>■ invisible), given by Eq. (11121) . In what follows, T(15') vv 
mode may be neglected in our analysis. 

Calculation of T(3S') — vv'y branching ratio within the Standard Model is straightfor- 
ward: it yields 

dB{T{3S) ^ VV^) _ N^Gl a /|(3S)^I(35) 

di " 16274^ 1^(35) ^ 

and 



24:671 47r i x (3?\ 



where Mt(3S) = 10.355^61^ and Press) = (20.32 ± 1.85)fce\^ y , The decay constant /t(3S) 
may be found, using 

P(T(35) - e+e-) = J^^''^ (05) 



where P(T(35) ^ e+e") = (0.443 ± 0.008)keV y . One gets /t(3s) = (0.430 ± 0.004)^61/ 
and, subsequently, 

B{T{3S) vvj) = (3.14l°:|) X 10"^ §iQ) 
12 



Thus, T{3S) — > z/z/7 branching ratio is about three orders of magnitude less than the BaBaR 
experimental limit on B{T{3S) ^ 7 + invisible), given by Eq. flllSp . In what follows, the 
effects related to T{3S) uu^y decay may be neglected as well. 

Thus, we may neglect the neutrino background effects when confronting theoretical pre- 
dictions for T decays into invisible states with the experimental data. 



4. MODEL-INDEPENDENT BOUNDS ON THE WILSON COEFFICIENTS 



In this section we use our theoretical predictions and existing experimental data to derive 
model-independent constraints on the Wilson coefficients Ci, 1=1,2,3,4, as functions of DM 
particle mass, m$, and the heavy mass. Ah- 

We start with the bounds, coming from the study of T(IS') $$* decay. These bounds 
will be used in the next section, to constrain the models with light complex scalar Dark 
Matter. As follows from Eq. f l2[5p . the decay branching ratio depends on the Wilson coef- 
ficient C3, the WIMP mass m$, the heavy mass Ah, as well as on the T(15') mass, total 
width, and decay constant. The numerical values of Mx(i5), rx(i5) and fr{is) have been 
specified in the previous section. Using these values, one gets 

3/2 



4 ^2 /lOOGeV\^ 



B{T{1S) = (5.3 ± 0.2) X 10-^ C| i^-^ j 



(HI) 



The uncertainty, about 4%, in the numerical factor in the r.h.s. of ( 14[ip is due to that in 
the values of the input parameters. Such a small uncertainty may safely be neglected during 
the further analysis. 

As it follows from Eqs. (131 2p and (I4[ip . T(IS') $$* branching ratio may be significantly 
greater than that of the T(IS') vv transition. In what follows, T(IS') decay into a 
pair of DM particles may be the dominant channel contributing to T(IS') invisible. 
Yet, in order to test this channel, the relevant experiments must be sensitive (at least) to 
Brr{lS) — > invisible) ~ 10~^. This sensitivity has been reached by the BaBaR experiment 
30|], as it follows from the bound on i?(T(15') invisible), given by Eq. f ll[2p . Substituting 
f ll[2p into fl4[ip . one derives the following constraint on IC3I as a function of m<j, and A^'- 
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FIG. 2: Upper bound on IC3I a) as a function of m$, for h.H = 100 GeV , b) as a function of Ajj, 
for = 1 GeV. 

The behavior of the upper bound on IC3I with the DM particle mass is presented in 
Fig. [2]^a). The bound on IC3I is almost insensitive to the WIMP mass for m$ < 2 GeV; 
it grows rater slowly for 2 GeV < m$ < 3 GeV; however as m$ > 3 GeV, the IC3I 
bound starts increasing rapidly, due to the rapidly shrinking phase space. For m$ < 3 GeV 
and Ah — 100 GeV, the derived bound on IC3I may be translated into constraints on the 
relevant couplings of models with a light complex spin-0 DM field. For m$ > 3 GeV, the 
experimental sensitivity to T(IS') — * invisible transition must be improved to compensate 
the phase space suppression. Experimental sensitivity improvement is also necessary for 
higher values of the heavy mass Ah'- the upper bound on IC3I grows rapidly with the heavy 
mass, as it can be seen from Fig. [2](b) and as can be inferred from the quadratic dependence 
of this bound on Ah- 

Note that within specific models, the value of G3 is somehow correlated to the values of 
the other Wilson coefficients, Gi, G2 and C4. In what follows, bound f l4[2p on IC3I may also 
lead to some constraints on Gi, G2 and C4, within particular models with a light complex 
spin-0 DM field. In terms of the branching ratios this means that experimental bound fll[2p 
on B{T{1S) invisible) (or on B{T{1S) $$*)) may lead to a phenomenological upper 
bound on B{T{3S) $$*7), within particular models with a light complex spin-0 DM field. 

At first glance, it may seem that a phenomenological upper bound on B{T{3S) —>■ $$*7) 
may have little practical use, in light of the existing BaBaR experimental constraint on 
B(T{3S) ■y+ invisible), given by Eq. f ll[3l) . Yet, constraint (11 [31) is derived for the emitted 
photon having monochromatic energy: it has been assumed that T(3S') — 7 + invisible 
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100 GeV 




FIG. 3: The differential branching ratio dB{T{3S) — > ^^'y)/ds versus the missing mass ^/s within 
a self-conjugate DM scenario for = 1 MeV (line 1), m<j, = 1 GeV (line 2), m$ = 2 GeV (line 
3) and m<j, = 3 GeV (line 4). 



transition is mediated by an intermediate resonant Higgs state A° 3l|]. Such a light Higgs 
state may exist e.g. within the extensions of the Minimal Supersymmetric Standard Model 
(MSSM) with an additional Higgs singlet Q, H, 39|. Bounds of ref. 31| on B{T{3S) ^ 
7 + invisible) have been plotted as a function of the mass of A^, or equivalently, of the final 
state fixed missing mass. 

In the case of non-resonant DM production considered here, when the decay is mediated 
by heavy degrees of freedom, the final state missing mass (or the photon energy at T rest 
frame) is not fixed. Instead, one has a broad missing mass distribution over the entire 
missing mass range, 2m$ < ^/s < Mx(3S)- For self-conjugate DM scenarios, the missing 
mass distribution shape is model-independent [as it is easy to see from Eq. (I2[2ip ] and is 
depicted in Fig. [3] for different choices of the WIMP mass. For complex scalar DM field 
scenarios, the missing mass distribution analysis depends on particular values of the Wilson 
coefficients [as one can see from Eqs. fl2[lip - fl2[13p ] and hence on a particular model. It 
is however clear that apart from the endpoints, it is non-negligible for the entire missing 
mass range as well. In other words, the experimental analysis, performed in 31||, should be 
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extended to the cases, when the emitted photon energy is non-monochromatic and is in the 
range < u < Mx(35)/2 — 2m|/Mx(3s)- To our knowledge, this work is in progress now^ 
Note that similar problems exist with the earlier bounds on T — 7 + invisible, reported 



by CLEO |33|. Bound 5(T(15) ^ 7 + X) < 3 x 10"^ is derived assuming that X is a 
single particle state. Thus, the emitted photon is monochromatic again. The only existing 
bound for the case of an invisible particle pair, B(T{1S) 7 + XX) < 10^'^, is too weak to 
produce any constraints within the models with non-resonant DM production. Due to the 
factor a/(47r) ^ 5.8 x 10""^, B{T $$*7) is always below the quoted limit. 

Of course, this all does not mean that the existing experimental constrains on B{T — >■ 
7 + invisible) are totally useless, if light spin-0 DM production is mediated by heavy non- 
resonant degrees of freedom. Note that BaBaR constraint f lllSp on B{T{3S) 'j+invisible), 
except for being plotted as a function of the final state missing mass, may also be rewritten 
as a bound for a missing mass interval. One can see from B{T{3S) — ^ 7 + invisible) versus 
mA_o plot of ref. jsil that 

B(T{3S) ^ 7 + invisible) < 3 x 10"^ (US) 

for y/s ^ 7 GeV or approximately s < and provided that the emitted photon 

energy is monochromatic. 

It has been discussed already that within the complex scalar DM field scenarios, one 
may derive a phenomenological upper bound on B{T{3S) —>■ $$*7). Comparison of this 
bound with fl4[3l) allows us to estimate if upcoming experimental constraints on T — 7 + 
invisible, for the case of non-monochromatic photon emission, may further improve the 
existing constraints on the parameter space of a considered model. We will use this approach 
in the next section. 

In this section, we concentrate hereafter on self-conjugate DM scenarios only. Recall 
that within these scenarios, T(IS') transition rate vanishes, thus we are left with 

T(3S') $$7 channel only. We perform the analysis of this decay channel, using the 
partially integrated branching ratio, for the missing mass interval, where bound fl4[3p is 
valid. We remind the reader that the insertion of a missing mass cut-off reduces significantly 
the WIMP mass range where an imposed experimental bound is efficient - we illustrate this 



^ The author is grateful to Yu. Kolomensky for this information. 



16 



Ci^ + C; 



16, Ah = 100 GeV 




2 3 
[GeV] 



FIG. 4: Partially integrated and total integrated branching ratios for T{3S) <I><I>7 decay (lines 
1 and 2 respectively), as functions of DM particle mass m<j,. The dashed line is the experimental 
bound 



in Fig. m Yet, as it has already been noted, all the existing bounds on T ^ 7 + invisible 
have been derived for a restricted missing mass range or for a restricted invisible particle 
mass range (much below the kinematical threshold). 



Using Sr^ 



x(35)/2 and [3J] Mr 



(35) 



10.355 GeV, Tr(3S) ~ 20.32 keV, /t(35) 



0.43 GeV, one may rewrite formula fl2[22p for T(3S') $$7 partially integrated branching 
ratio in the following form: 



B(T(3S) ^ **7)|«M»,„,/2 



2.6xlO-'(Cr + C|)(12^)VM 



H 



where 



/(x$) 



4 



1 + ) y/T 



32 



— — x$(l 



x<s,) In 



1 + VI - 8x$ 
2^2. 



(!14) 



(!15) 



Note that < /(x$) < 1, /(x$) = 1 if x$ = 0, and /(x$) = if x<i, = ml/M^^-g^-, 



1/8. 



At first glance, it may seem from Eq. f l4[4l) that T{3S) $$7 branching ratio is far out 
of reach of the BaBaR experimental sensitivity, for a reasonable choice of Ci and C2. Notice, 
however, that within certain models with light spin-0 Dark Matter, the Wilson coefficients 
Gi and/or G2 may be enormously large, as they contain some enhancement factors, such 
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FIG. 5: Upper bound on JCf + Cl as a function of m$, for Ah = 100 GeV. 



as the ratio Ku/^b ^ 1 - due to the mass term in the numerator of a heavy fermion 
propagator, or the Higgs vev's ratio tan/5 (with the latter being, say, ~ mj/mfe ^ 1) - due 
to DM particle pair production via exchange of a heavy non-SM Higgs degree of freedom. 

These enhancement factors can make Ci and/or C2 to be ^ 10 and hence -B(T(3S') — *• 
$$7) to be ~ 10^^ — 10~^, i.e. significantly exceeding bound (14131) on B{T '-y + invisible). 
Yet, bound (I4[3I) assumes that the emitted photon energy is monochromatic: rigorously 
speaking, one should wait until the experimental limit on B{T — > '-y + invisible) for the case 
of non-monochromatic photon emission comes out^. One may use (I4[3I) only to derive a pre- 



liminary estimate of possible constraints on y Ci + and hence on the relevant parameters 
of light spin-0 self-conjugate DM models. 
This estimate may be presented as 



^ Also, the limit may be derived for T(IS') state instead of T(3S') - the author thanks Yu. Kolomensky 
for the discussion of this point. The reader, however, can easily check that making the replacements 
Afx(3S) My(is)i /t(3S) ^ It (IS) ^'iid rT(3S) ^ rT(is) formula ()2I22|) modifies our predictions by 
about 25% only. 
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a) b) 

FIG. 6: Tree level diagrams for T(15) decay within New Physics models. The transition 

occurs a) by exchange a gauge boson, b) by exchange of a charge -1/3 fermion. 

It is also depicted in Fig. [5], as a function of m$, for Ah = 100 GeV. Based on the discussion 
above, this estimate implies rigorous constraints on the relevant parameters of the models 
with light spin-0 self-conjugate Dark Matter, for m$ < 3 GeV . We analyze these constraints 
within a particular model in Section [61 



5. COMPLEX DM FIELD SCENARIOS: MIRROR FERMION MODELS 

As it was mentioned above, the distinct feature of the scenarios with a light complex 
spin-0 DM field is that T(IS') $$* decay rate is non- vanishing. One may use within 
these scenarios bound fl4[2p on IC3I, both to constrain the model parameter space, and - due 
to possible correlations in the values of the Wilson coefficients - to derive a phenomenological 
upper bound on the T{3S) $$*7 branching ratio. 

Recall that bound (141 2p on IC3I is strongest for Ah — 100 GeV and m$ < 3 GeV. Yet, 
even for these values of the heavy and WIMP masses, C3 is still allowed to be of the order 
of unity. It seems to be very unlikely to saturate such a (rather weak) bound, if within a 
full electroweak theory, T(IS') $$* transition is loop-induced. We may therefore restrict 
ourselves by the models, where T(IS') $$* decay occurs at tree level. 

Possible tree level diagrams for T(IS') $$* transition within New Physics models are 
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presented in Fig. [61 They involve exchange of either a neutral gauge boson or a charge -1/3 
fermion. These type of diagrams may occur e.g. within U-boson models^ or within mirror 



fermion models with light scalar Dark Matter [6|, 



MM- 



Models with a new (beyond the SM) neutral gauge U-boson have been considered yet 



ong ago, both within the supersymmetric theories and within the other SM extensions 



40| 



44( 1 . T decays into invisible states (with or without a photon emission) have been studied 



in details in 







I, |45|, within the scenarios with a light U-boson (with a mass less than a 
hundred MeV scale). In this paper only scenarios with a heavy U-boson are of interest. In 
that case, however, U-boson should be much heavier than the Z Both T(IS') $$* and 
T(3S') $$*7 branching ratios (inversely proportional to mfj) would be then suppressed 
too much to yield any constraints on the parameters of U-boson models. 

Thus, we concentrate here only on the models with mirror fermions, where T meson 
decays into Dark Matter by exchange of a heavy, charge -1/3 fermion, as shown in Fig. [6](b). 
Within these models, scalar Dark Matter couples to ordinary matter by means of Yukawa 

tin 

interactions, such as p, 113[ 

- £ = $ [Xb^ Fb^ bL + Afe^ Fb^ Br) + h.c. + ... Pl) 

where Fh is a charge -1/3 colored mirror fermion. Unlike ordinary quarks, however, the right- 
handed component of Fb transforms as a member of a weak isospin gauge doublet, while 
the left-handed component of Fb transforms as a singlet. In other words, Fb appears to be 
a mirror counterpart of b-quark. Fb and other mirror fermions do not mix with ordinary 
quarks and leptons, as they - along with the DM particle - are odd under so-called M-parity 
Gj, or mirror parity transformation, whereas ordinary matter is M-even. The scalar DM 
particle is the lightest M-odd particle of the theory. 

Mirror fermion scenario may be realized [l3, Q, [l^ within the MSSM with gauge- 
mediated supersymmetry breaking and the DM particle in the hidden sector. It has been 



argued 13| that the thermal relic density constraint may be satisfied irrespectively of the DM 
particle mass and, in particular, within the light DM scenario (the "WIMPless miracle"). 
As for the mirror fermions, they serve as connectors between the hidden and the MSSM 
sectors. For the sake of simplicity, it is assumed that each connector couples to one quark 



^ Scenarios with T decaying into DM with a large rate through the SM Z boson may be excluded right 
away, as it would imply that Dark Matter contributes to the Z boson invisible width. 
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(or lepton), but other scenarios are in principle possible, where each mirror fermion has 
multiple couplings with the SM fermions, or each quark or lepton couples to multiple mirror 
fermions. The masses of mirror fermions are expected to be of the order of electroweak 
braking scale or heavier^, i.e. Mp ^ 100 GeV . 

The tree level matching between the full and effective theories yields 



Ci = - 



+ 



bLl 



|2 



ml J 2 4 



Co = - 



ml 



^ _ \XiJ^ + \>^ ^ A^ = M^, (02) 



where m^ is the MS bottom mass evaluated at the matching scale. The natural choice 
of the matching scale is the heavy mass, fin = ^h, or, accor ding (151 2p . fin = Mp^. The 



bottom mass evolution with the scale is known up to four loops j48j, yet to the leading-order 
approximation used here, it is given by 

ml = m,{MpJ=m,{m,)(^^!^^] ' (03) 



where m^,{m^) = (4.20 ± 0.1 7) Gel/ [34]. 

Running of with the scale is the only 0(1) QCD effect, relevant for our analysis. As it 
has been already mentioned above, the Wilson coefficients Gj are low-energy renormalization 
scale-independent. That is to say, equations (151 2p hold also at any scale fi < Mp^, including 
the decay scale. 

We begin with using (151 2p to transform bound (14120 on IG3I, as a function of m$ and 
Ah, into that on the couplings Aj,^ and Af,^, as functions of m$ and Mp^. This bound, in 
general, depends on possible correlations in the values of A^^ and Af,^. Here the following 
two scenarios are considered: 

• the chiral scenario, when one of these coupling vanishes, e.g. A^^ = 0; 

• non-chiral scenario, with A;,^ = A;,^ = A;,. 



Unlike 

we do not apply the bound on the fourth generation quark mass, nidi > 258 GeV 
47[. In our opinion, this bound is irrelevant for the mirror fermion models. Indeed, constraints coming 
from the annihilation channel, d^d^ qqWW, are invalid here: mirror fermions may annihilate to a DM 
particles pair (through interactions like that in (|5lip ) and hence escape detection. The other channel in 
use, ^4 cW, is also invalid, due to the mirror symmetry of the model. 
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Within the chiral scenario, the use of fl4[2l) and (151 2p yields 

To the best of our knowledge, bound on the coupling A;,^ (equivalently on A^^, if we choose 
instead A;,^ = 0) is derived for the first time. Dark Matter direct search experiments, based 
on DM scattering off nuclei, are sensitive to the light quark - mirror fermion interaction 
couplings only 13|, [iJ] . So is the dominant contribution to DM annihilation processes. The 
decays B K + invisible or Bg — * invisible would inevitably depend both on A^^ ^ and 
on As^^, but not on A;,^ ^ alone. Other quarkonium, Xbo, invisible decay modes still need 
improvement of experimental sensitivity 46|. 

At first glance, bound fl5[4p on A;,^ is weak. Furthermore, it is essential only for a restricted 
range of the mirror fermion mass: it becomes weaker than the perturbativity limit, Xb^ < 
\/47r, if Mp^ ^ 200 GeV. Nevertheless, the use of f l5[4p may lead to phenomenological upper 
bounds on some of bottomonium decay channels, and in particular, on T{3S) $$*7. 

Indeed, for the chiral scenario, matching conditions fl5[2p may be rewritten in the following 
form: 

Ci=C3 = -C, = C2 = (05) 

Thus, bound f l5l4p on | A;,^ | may be transformed into the constraints on the Wilson coefficients 
Ci and C4. Then, using formulae f l2[18p - (121 20p . for C2 = and the other parameters values 
specified in Sections 2 and 3, and applying the constraints on Ci and C4, one derives an 
upper bound on T(3S') $$*7 (total) branching ratio, as a function of DM particle mass, 

We present this bound in Fig. [71 It may also be rewritten as -B(T(3S') $$*7) < 
1.57 X 10~^, with the limit being saturated, as m$ 0. Certainly, within the chiral 
scenario, B(T{?)S) $$*7) is far away of the present experimental sensitivity. This result 
is not surprising: it has been anticipated in the previous section, for the models with no 
enhancement factors in the Wilson coefficients. The use of the derived constraint fl5[4p on 
lAft^l enables one to formulate this quantitatively, by imposing a well defined limit on the 
T(3S') $$*7 branching ratio. 

It is worth noting that a signal for the T(3S') ^ 7 + invisible transition, above the 
quoted limit on i?(T(3S') $$*7), would rule out the mirror fermion models with light 
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FIG. 7: Phenomenological upper bound on B(T{3S) ^^*^) branching ratio within the chiral 
scenario. 

complex spin-0 Dark Matter and chiral mirror couplings. The following statement is also 
true: such a signal would imply that light spin-0 DM field is self-conjugate within the mirror 
fermion models with chiral mirror couplings. This result is not surprising as well: it has 
been noted previously that study of T ^ $$* transition (combined with that of T — $$*7) 
represents an opportunity to test, whether Dark Matter, if being light and a spin-0 state, is 
self-conjugate or complex. 

Within the non-chiral scenario, A^^ = A;,^ = Af,, noting that C3 = — |Abp/2, we may 
rewrite bound (14121) on IC3I in the following form: 

(\ -3/8 

We present this bound in Fig. [HI as a function of m$, for Mp^ = 100 GeV. The derived 
constraint on |Afe| changes rather slowly with the WIMP mass, if m$ < 3 GeV. To simplify 
the analysis, one may use, within a crude approximation, \Xh\ < 1.22 for this range of the 
WIMP mass and Mp^ — 100 GeV. The behavior of bound (15161) with the mirror fermion 
mass is rather trivial (max[|Afc|] grows linearly with MpJ. We note here only that (I5[6p 
should be replaced by the perturbativity limit. A;, < \/47r, if Mp^ > 300 GeV. 
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FIG. 8: Upper bound on |Ab| as a function of m$, for Mp^ = 100 GeV. 



Generally speaking, it is expected that ordinary-to-mirror fermion couplings be signif- 
icantly le. than one fl Q. It the=e couplings a.e ot the san.e orie. for all th.ee quark 

generations, one expects A;, ~ A„,d ~ 0.1 (Mp^^/lOOGeV) [6| from the DM relic abundance 
condition, Q^Mh"^ ~ 0.11. Otherwise, if a hierarchy in the couplings exists, one may argue 
that this expectation does not hold for A;,: its value is not affected by the Dark Matter direct 
search designated experiments, nor does it have essential impact on DM annihilation and, 
hence, relic abundance. Furthermore, having a third generation Yukawa coupling to be of 
order unity is not unusual. However, bound (151 6p on |Ab| may be significantly improved due 
to possible constraints coming from the study of T ^ $$*7 transition. 

Indeed, within the non-chiral scenario, matching conditions (151 2p may be rewritten for 
the Wilson coefficients Ci, C2 and C4 as 



Co — Ca 







P7) 



Because of the enhancement factor Mp^/ml, Ci may be enormously large: as it has been 
discussed in Section HI T $$*7 branching ratio is then within the reach of the present ex- 
perimental sensitivity. Furthermore, due to this enhancement factor, B{T{3S) —>■ $$*7) oc 
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roughly^, unlike B{T{1S) $$*) oc Thus, within the non-chiral scenario, 

T(3S') $$*7 mode may be sensitive to a wider range of the mirror fermion mass than 
T{1S) 

For the numerical analysis, we use lA^I < 1.22. As it is mentioned above, this bound may 
be used in a crude approximation for m$ < 3 GeV and Mp^ = 100 GeV. In order to highlight 
the behavior of B{T{3S) $$*7) with the mirror fermion mass properly, we keep using 
|A;,| < 1.22 for Mp^ > 100 GeV as well (even though the constraint on |Afe| is significantly 
weaker then). This yields 5 (T (35) ^ $$*7) < 1.23xl0-^ B{T{3S) <l'$*7) < 3.4 xlO-^ 
and B{T{3S) <l><l>*7) < 0.93 x IQ-^, for Mp^ = 100 GeV, Mp, = 200 GeV and Mp, = 
400 GeV respectively, with the limits being saturated, as m$ — ^ 0. We compare these 
limits to bound fl4[3l) on i?(T(35') 7 + invisible). Recall that bound fl4[3p is derived for 
the emitted photon having monochromatic energy, thus it may be invalid for the scenarios 
where light DM is produced due to exchange of heavy non- resonant degrees of freedom. Yet, 
comparison of the above derived limits on i?(T(35') —>■ <l>$*7) to fl4[3p allows us to infer 
that upcoming BaBaR constraint on T — ^> 7 + invisible, for the case of non-monochromatic 
photon emission, may essentially improve the existing bound on |A;,|. 

If this constraint is of the same order as fl4[3p . one would get lA^I ^ 0.5, |A;,| ^ 0.65 
and \Xb\ < 0.9, for sufficiently low WIMP mass and for Mp^ = 100 GeV, Mp^ = 200 GeV 
and Mp^ = 400 GeV respectively. The upcoming experimental data, however, may be not 
so optimistic, they may lead to a much weaker experimental limit than fl4[3p . In either 
case, we want to emphasize again that we derive here preliminary estimates only of possible 
improvements of the existing bound on |A;,|. More rigorously, regardless of any expectation, 
one must presently use the existing bound on |Aft|, given by Eq. (151 6p . 

To summarize our discussion of the mirror fermion models, we note that study of T(IS') 
decay into a spin-0 DM particles pair, T(IS') $$*, leads - for the first time - to the 
bounds on the couplings of b-quark interaction with its mirror counterpart. Within the 
non-chiral scenario, these bounds may be somewhat improved by upcoming constraints on 
the T — i> 7 + invisible mode for a non-monochromatic photon emission. Recall also that we 
assumed throughout this section that light spin-0 DM field has complex nature. Otherwise, 
within the self-conjugate DM scenario, no constraints on A^ couplings, apart from preliminary 

^ Slight deviation from this rule occurs due to running of the bottom mass with the heavy mass scale. 
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estimates, may presently be derived. 



6. SELF-CONJUGATE DM SCENARIO: DARK MATTER MODEL WITH TWO 
HIGGS DOUBLETS (2HDM) 



In this section we will discuss the models, where interaction of (self-conjugate) spin-0 
Dark Matter with the ordinary matter is mediated by heavy Higgs degrees of freedom. 
These models are known to be the most economical SM extensions, containing Dark Matter: 
they are created by extending (if necessary) the SM Higgs sector and embedding a scalar 
DM particle into theory, by adding a rather small number of unknown parameters. 



The simplest model of this type is the Minimal Scalar Dark Matter Model |10l. I49l. |50[|. It 
has the same particle content as the SM, plus a gauge singlet real scalar field odd under 
discrete symmetry ($ — $) and coupled to the SM particles through the exchange of Higgs 

-iSl^ Presently, 
53|, if the 
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boson. This model has been widely studied in the literature |lOl . 
the Minimal Scalar DM Model has a very restrained parameter space 
DM particle is chosen to have a GeV or smaller mass. Also, the T decay channels considered 
here are not sensitive to the parameter space of this model. Indeed, as it has been shown 
in Section HI for i?(T(3S') $$7) to be within the reach of the present experimental 
sensitivity, the Wilson coefficients C\ and must contain some enhancement factors, such 
as heavy-to-light mass ratio or a large Higgs vev's ratio. None of these enhancement factors 
may be generated in the Minimal Scalar Dark Matter Model, where DM particle production 
occurs solely via the SM Higgs boson exchange. Recall also that the other decay, T — > $$*, 
rate vanishes, when $ = $*. Instead, the Minimal Scalar DM Model is well tested by 
B ^ K + invisible mode, if assuming that DM particle has a GeV or smaller mass. Study 
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12| 



of this mode leads to rigorous bounds on the model parameter space 

In this section we consider the simplest extension of the Minimal Scalar DM Model, the 
two-Higgs doublet model (2HDM) with a gauge singlet real scalar DM particle. The DM 
interaction part of Lagrangian, relevant for our analysis, may be written as jll] 



Tfl 



(01) 
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where $ is a odd real scalar DM field and 




Ho 



2 I U U IjO rrO tt+ tj— 



Following Refs. |62| , we consider type II version of 2HDM, where Hi generates masses 
of down- type quarks and charged leptons, whereas H2 generates masses of up-type quarks^. 
The Higgs vev's, fi, f2, are constrained by the condition v1 + = v"^ = {2AQGeVY . The 
Higgs vev's ratio, tan/3 = ^2/^1, is a free parameter of the theory. We assume here that 
tan/? ^ 1. As discussed before, i?(T(3S') $$7) is then enhanced and may be within 
the reach of the BaBaR experimental sensitivity. Large tan (3 scenarios are also known to 
be much less restrained than the Minimal Scalar DM Model. One may avoid tension with 
satisfying the DM relic density constraint. Another point worth mentioning is that the 
fine-tuning of parameters needed to generate a GeV or smaller DM mass is, in general, 



significantly weaker llj than that within the Minimal Scalar DM Model. 



After eliminating Goldstone modes, one may express Hi and H2 in terms of two CP-even, 



one CP-odd and two complex charged mass eigenstates - h^, H^, H"^ respectively |63| . 



6^: 



H^ = ^(vi + H° cosC - /i° sine + sin /?) , H+ = H+ sin /3 



V2 

H^ = ^[v2 + h°cosC + H°smC + iA^cosp), H^ = H- cos P P2) 

Here and H^ are the lightest and heaviest CP-even states respectively, and ^ is the 
CP-even mass-matrix diagonalization angle. 

Using Eqs. (161 2p . one may rewrite (I6[ll) in the following form: 

2 

-C = + A;,o$ V $2/1° + X^,^ V $2/7° + . . . P3) 

where 

,2 



+ (Ai cos^ /? + A2 sin^ 13 + 2A3 sin (3 cos (3) (04) 
is the DM particle mass, 

•^ho$ = — Ai sin^ cos /? + A2 cos^ sin /3 + A3 cos(e + /?), (jHlS) 

-^_ffO# = -^1 cos^ cos /5 + A2 sin^ sin /? + A3 sin(e -|- /3), ([6l6) 



^ The used definition of Hi and II2 corresponds to the following notations for the Yukawa interactions: 



27 



and the ellipsis in fl6[3p stands for the dropped quartic interaction terms. To the leading 
order in the perturbation theory, T(35') $$7 transition occurs at tree level by exchange 
of a single Higgs boson. Thus, to the leading-order approximation, it is sufficient to consider 
only the cubic interaction terms in fl6[3p . Then, as it follows from fl6[3p . only the CP-even 
Higgs states are relevant for our analysis. 



The b-quark Yukawa interaction terms may be written in the following form 3J, |62|, |63| : 



_Cy = -^^-bh'b + ^'-^bH'b + ... (07) 

V COS p V COS p 

where we write down explicitly only the b-quark interactions with the CP-even Higgs bosons. 
Our strategy is the following now. We use fl6[3p . (161 5p - fl6[7p . to derive the matching con- 



ditions for the Wilson coefficients Ci and C2. Then, we transform bound fl4[6p on yCf + C| 
into that on the relevant couplings of the model, depending on the DM particle mass, Higgs 
mass and tan (3. We have to recall, however, that this bound may serve only as a preliminary 
estimate of possible constraints on the parameter space of the model. 

In the limit of tan/3 ^ 1, the CP-even mixing angle ^ has two possible solutions 

(a) e ~ vr/2 - /3 

(b) i ^ -/5 

In the former ), the lightest CP-even Higgs boson, , is Standard Model-like: 

its phenomenology is similar to that of the SM Higgs boson and the experimental bound on 
its mass is close to the SM limit'' (see [3^ and references therein). 

In the latter case (case b), is "New-Physics (NP) like": its phenomenology differs 
drastically from that of the SM Higgs boson [63|. In particular, m^o may be much below 
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67|, 



the Standard Model experimental limit: according the existing experimental data 
m^o > 55 GeV or 171^0 < 1 GeV in the general type II 2HDM. 

As it was mentioned above, the light Higgs scenario is beyond the scope of the present 
paper, thus we assume here that rriho > 55 GeV. Notice, however, that this bound is derived, 
provided that no invisible Higgs decay mode exists. On the other hand, if the NP-like Higgs 
invisible decay mode is dominant, it may escape detection. No bound on m/jO, to our best 
knowledge, exists in that case. 



Also, if the SM Higgs decays predominantly invisibly, the SM lower experimental bound is distorted by a 



few GeV only |M,l65|. 
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Within the considered model with hght scalar Dark Matter, analysis of the T $$7 
mode may restrict the scenarios with an invisibly decaying lightest Higgs boson by putting 
severe constraints on the interaction coupling, A/jO^,, given by Eq. fl6[5p . If the lightest 

CP-even Higgs boson is NP-like ~ —jS, case b), one may also rewrite fl6[5l) as 

A/^o$^ A3 + (Ai + A2)cos/? P8) 

The last term in the r.h.s. of (161 Sp . although being suppressed by a factor of cos /? ~ 1/ tan /3, 
must be retained because of possible hierarchy in the values of A3 and Ai or A2. Scenarios 
with such a hierarchy may be of importance, as A/tO$ is constrained to be 0(1/ tan/3), if /i° 
mass approaches to its lower limit, 171^0 = 55 GeV. 

Indeed, neglecting the heaviest SM-like CP-even Higgs exchange contribution, the match- 
ing conditions for the Wilson coefficients have the following form: 

Ci = tan/3, C2 = 0, Ah = rriho (09) 

Thus, one may rewrite bound (I4[6p on JCf + as 




(010) 



where f{x^), Xq> = rn%/M'^^^gy is given by Eq. (I4[5I) . 

Formula (I6[10p implies a stringent upper bound on \Xh»q>\ at the Higgs mass lower thresh- 
old, ruho = 55 GeV. Choosing e.g. tan P = 20, one has |A/iO$| < 0.1 for the WIMP mass, m$, 
less than 2 GeV, as one can see from Fig. [91 For the other choice of the Higgs vev's ratio, 
tan/3 = 40, the same constraint on \XhO^\ is derived for m$ = (2.5 — 3)GeV, one also gets 
|AhO$| < 0.05, as m$ < 1 GeV. For such small values of |A/jO$|, it seems to be very unlikely 
that would escape detection and its mass be below 55 GeV. More rigorously, however, 
detailed reanalysis of the Higgs production and decay rates, including that of —>■ $<!', 
should be performed (which is beyond the scope of the present paper). 

Bound on |A/iO$| may still be rigorous, if the mass is heavier than 55 GeV. For instance, 
choosing m/jo = 100 GeV (in which case hf^ is still NP-like) and m$ < 2 GeV, one gets 
|AftO$| < 0.2 for tan/3 = 40, m$ < 1.5 GeV and for tan/3 = mt{mt) / mf,{mt) ~ 55, m$ < 
2 GeV. Yet, Xho<s> constraints become weaker with decreasing tan/3 and increasing WIMP 
mass. It may be of order of the SM weak coupling, if tan/3 ~ 20 and m$ ~ 2.4 GeV. 
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FIG. 9: Upper bound on |A/jO$| as a function of m$, for mf^o = 55 GeV and tan/3 = 20 (line 1), 
tan/? = 40 (line 2). 

Thus, within the type II 2HDM with a light spin-0 Dark Matter, study of T — $$7 
decay channel may lead to severe constraints on the lightest CP-even Higgs invisible decay 
coupling, if that Higgs is New-Physics hke. 

As it follows from Eq. fl6[8l) . bound fl6[10p on |A/iO$| may also be transformed into that on 
the couplings Ai, A2 and A3. Constraints on Ai and/or A2 may also be derived from the study 



oi B ^ K + invisible transition 



Those, in general, are strong enough: in particular. 



0(1) values of Ai and A2 for niho ~ 55 GeV are ruled out, if m$ < 1.5 GeV. In that case, 
1/ tan/3 suppressed terms in Eq. (16 [8 1) may safely be neglected, so that A/iO$ ^ A3. In other 
words, bound (I6[10p on |A/iO$| is also that on IA3I, if m$ ^1.5 GeV. 

Note that due to cancellation effects in the relevant diagrams, WIMP pair production 
rate in B meson decays is insensitive to the value of A3 [11], and hence, in general, to the 
value of Xho<i>. The scenarios with A3 dominant, or at least non- negligible, have been thus 
far unconstrained. To our best knowledge, bound on the interaction coupling, A/tO<j,, is 

derived for the first time. Thus, study of DM production in T decays enables one to test the 
regions of the parameter space of 2HDM with scalar Dark Matter, which are inaccessible by 
B meson decays with missing energy. 
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It may seem naively that the rigorous constraint on the interaction couphng, Xho<s,, 

results also to suppress the DM annihilation rate, which in its turn may lead to scenarios 
with overabundant Dark Matter. Yet, the DM annihilation rate to down-type quarks and 
charged leptons is enhanced by a factor of tan^ (3 (due to tan (3 enhancement of these fermions 
Yukawa interaction with the CP-even Higgs, if the one is NP-like). Or, equivalently, using 



the notations of ref. 



the effective coupling for these annihilation processes is 



flOOGeV\^ 

Kyfiq> = XhOij, tan p 

V rriho / 

In terms of the coupling k^o^, bound ( I6[10p may be rewritten as 

Kf,o^ < 6.8 r'/\x^) Pll) 

Keeping in mind that < f{x^) < 1 or, equivalently, 1 < /^^^^(a;$) < oo, bound (161 lip on 
K,ho<i> yields no essential constraints on DM annihilation and, hence, no tension with satisfying 
.DM relic abundance constraint above the Lee- Weinberg limit of the model, m$ ^ 100 MeV 



Thus, the 2HDM scenario with A3 dominant, or at least non- negligible, is much less 
restrained, than those with Ai or A2 dominant or than Minimal Scalar DM Model jll. [l^. 

So far it has been assumed that the lightest CP-even Higgs boson is NP-like. In the 
opposite case, when is the SM-like (case a), the matching between the full and effective 
theories yields 

^ n '^i 1 / A2 A3tan/3\ 

In deriving (161 121) . we used, for tan/5 ^ 1 and ^ ~ 7i/2 — /3, 

Aho# A2, A/^o$ ^ A3 (El 13) 

It is not hard to see from our further analysis that omitted 0(1/ tan/5) terms in Eqs. (161 13p 
are not essential in this case. 

The first term in the expression for Ci/K\j is due to the SM-like Higgs exchange. As 
expected, it has no enhancement factor - thus, as it was discussed in Section HJ contribution 
of this term to the T(3S') $$7 rate is by (at least) an order of magnitude lower than 
the present experimental sensitivity. We may further disregard the Dark Matter interaction 
with the lightest CP-even Higgs h^. Then, we may rewrite matching conditions (161120 in a 
more transparent form: 

Ci = ^^^^, C2 = 0, AH = mHO P14) 
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FIG. 10: Upper bound on | A3 1, a) as a function of m$ for mjjo = 160 GeV and tan/3 = 30 (line 
1), tan/3 = 55 (line 2), b) as a function of m^jo for tan/3 = 55 and m$ = 100 MeV (solid line), 
m$ = 1.5 GeV (dashed-dotted line), m$ = 2 GeV (dashed line). 

In other words, we restrict ourselves by considering the contribution to T(3S') $$7 
amplitude due to exchange of the heaviest (NP-like) Higgs boson only. This contribution 
is enhanced by tan/5 factor, coming from hH^h Yukawa interaction coupling in (161 7p (if 
^ — (3). The remarkable feature of the considered scenario is that even though 

T(3S') $$7 transition is generated by exchange of the heaviest CP-even Higgs boson, 
the decay branching ratio for a certain range of if" mass is well within the reach of the 
present experimental sensitivity, due to tan^/3 enhancement. As a consequence, one may 
derive constraints on the coupling A3, as a function of m^^o, m<^ and tan/3. 



Indeed, using Eqs. (161 14p . one may rewrite bound fl4[6p on JCl + Cl as 



IA3I < 



17.4 



(- 



771^0 



dSis) 



tanpj \160GeVj 

Choosing mj^o = 160 GeV as a reference value is not accidental. Within the general type II 



2HDM, theoretical upper bound restricts the SM-like Higgs to be less than 180 GeV 68|. 
Also, the SM Higgs mass interval (160 — 170) GeV has been recently excluded with 95% C. 



L. by the CDF and DO data 



69j. Thus, within type II 2HDM, above 160 GeV, the CP-even 



Higgs boson is presumably the heaviest one and NP-like. 

As one can see from Eq. (161150 and Fig. [TOl^ a). for rriHO = 160 GeV and tan/3 = 30, A3 is 
constrained to be of order of the SM weak coupling or smaller (IA3I ^ 0.65), if m$ ^ 1 GeV. 
Also, for the same choice of the Higgs mass and tan/3, IA3I is to be less than one, if the 
WIMP mass is less than 2 GeV. Bound on | A3 1 is significantly more rigorous for higher values 
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of tan/3. For instance, if choosing tan/3 = mt{mt) / rribimt) ~ 55, one gets IA3I ^ 0.35 and 
IA3I < 0.5 for m$ < 1 GeV and m<j, ^ 2 Gel^ respectively. 

The restrictions on IA3I are essential also for higher values of the heaviest CP-even Higgs 
mass: for tan/5 = 55, they are still of interest up to m^^o ~ 280 GeV, as one can see from 
Fig. [TUT b). Our analysis may be spread for the uiho < 160 GeV range as well, leading to 
more rigorous constraints than those in Fig. [TOl 

Recall that bound on IA3I is derived for the first time. It is also worth noting that 
constraints on IA3I imply also those on the heaviest CP-even Higgs invisible decay rate, if 
that Higgs is NP-like. As it follows from Eq. (161131) . A3 is the interaction coupling in 

this case. Yet, as in the case (b), bound fl6[15p does not seem to have an impact on the DM 
annihilation rate (which is tan^/3 enhanced). Vice versa, DM annihilation processes (and 
possibly scattering off nucleons) may lead to additional constraints on A3. Study of these 
processes, however, goes beyond the scope of the present paper. 

Thus, within type II 2HDM with a scalar Dark Matter, for large tan (3 scenario, T meson 
decay into a Dark Matter particles pair and a photon, T $$7, may be used to derive 
essential constraints on the parameters of the model, which otherwise cannot be tested by 
B meson decays with invisible outcoming particles. 

7. CONCLUSIONS AND SUMMARY 

Thus, spin-0 Dark Matter production in T meson decays has been investigated. We 
restricted ourselves by consideration of the models where the decays occur due to exchange 
of heavy non-resonant degrees of freedom. Both the scenarios with a complex scalar DM 
field and those with DM particle being its own antiparticle have been analyzed. 

We performed our calculations within low-energy effective theory, integrating out heavy 
degrees of freedom. This way we derived model-independent formulae for the considered 
branching ratios. We used these formulae to confront our theoretical predictions with exist- 
ing experimental data on invisible T decays, both in a model-independent way and within 
particular models. It has been shown that within the considered class of models, DM pro- 
duction rate in T decays is within the reach of the present experimental sensitivity. Thus, 
T meson decays into Dark Matter, with or without a photon emission, may be used to con- 
strain the models with a GeV or lighter spin-0 DM. In particular, within the mirror fermion 
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models, using the existing BaBaR constraint on the T(15') — > invisible mode, we derived 
for the first time bounds on the parameters of the model that otherwise could not be tested 
by other DM search processes. 

Experimental constraints on the other mode, T{3S) — > 7 + invisible, are derived assum- 
ing that Dark Matter is produced by exchange of a light resonant scalar state. Within the 
scenarios with non-resonant DM production, these constraints may be used only to make 
preliminary estimates of possible bounds on the parameters of the models. Yet, those es- 
timates show that these bounds may be rigorous enough; besides, they are derived within 
the least restrained presently light scalar DM scenarios. Our goal is thus to encourage the 
experimental groups to analyze the experimental data on T ^ 7 + invisible also for the 
case of non-monochromatic photon emission and spin-0 invisible states. 

So, from our analysis one may conclude that Dark Matter production in T meson decays 
may serve as an interesting alternative to commonly used DM search methods, capable of 
providing a valuable information on DM particles, if those turn to have a mass of the order 
of a few GeV or smaller. 
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